Interference effects in resonant magneto-transport 
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We study non-equilibrium magneto-transport through a single electron transistor or an impurity. 
We find that due to spin-flip transitions, generated by the spin-orbit interaction, the spectral density 
of the tunneling current fluctuations develops a distinct peak at the frequency of Zeeman splitting. 
This mechanism explains modulation in the tunneling current at the Larmor frequency observed in 
scanning tunneling microscope (STM) experiments and can be utilized as a detector for single spin 
measurement. 
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Problems involving transport through low dimensional 
structures have received significant attention in connec- 
tion with the rapid development of spintronic 0] and sin- 
gle electron |2) devices. Apart from being likely candi- 
dates for becoming components of future electronic in- 
tegrated circuits, the use of single electron transistors 
(SET) as charge detectors have been contemplated in 
several solid state quantum computing designs 

In the present work we demonstrate that the spectral 
density of current fluctuations of a single electron tran- 
sistor in the external magnetic field develops a peak at 
the electron Zeeman frequency generated by spin-orbit 
interactions. We attribute such effect to the interference 
between the spin up and spin down components of the 
transmitted current resulting from the spin flips in the 
tunneling process. 

As a model system we consider a heterostructure (for 
example Si/Ge) schematically shown in Fig. 1. The two 
regions, to the right and to the left from the dotted line 
denoting the interface, have different ^-factors, « 2 for 
the left region and (72 7^ 2 for the right region. There are 
two contacts/Fermi reservoirs in each of the regions. The 
left region also contains a quantum dot, so that when a 
potential difference V is applied between the two reser- 
voirs, electrons can tunnel from left to the right reservoirs 
via the dot. The energy levels of the dot are spin-split 
by an external magnetic field. In this case the spin-orbit 
coupling causes the spin-flip transitions resulting in co- 
herent effects in the tunneling current jdj . 

We describe our system by the Hamiltonian H = 
Hl + Hji + Hs + He + Ht where the first two terms 
represent the unperturbed states of two contacts, = 
Y.i,s^isa\^°^is and Hr = Y.r.s'^^^^'lsO'rs, where (4^) 
creates a fermion/electron at the energy level e/ (e^) and 
with spin s in the left (right) reservoir. We assume that 
there is a single discrete level in the dot due to spatial 
quantization. The level is spin-split by the magnetic 
field B, so that the states in the dot are described by 
Hs ~ £s^s, where fig = aja^, and a\ creates an elec- 
tron in the dot at the level with spin s. We denote 
e_i/2 ei/2 = gPS = E, Fig. 1, where g is the electronic 



(^-factor in the dot and [3 is Bohr's magneton. The term 



corresponds to the Coulomb charging 



He = Es 

energy for the electron in the well. In what follows we 
will assume the case of complete Coulomb blockade, i.e., 
U 00, thus allowing for only one electron to occupy 
the two spin states in the dot. 
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Fig. 1: Quantum dot coupled to two contacts. The 
right contact has g-factor different from that of the left 
contact and the dot. Tunneling with spin flip generates 
effective coupling between the two Zeeman sublevels in 
the dot. 

The tunneling transitions between the left reservoir, 
the dot and the right reservoir are represented by the 
Hamiltonian: 



Ht ^"^^i (aL"s + alais 



+ ^ rirss' [almas' +al,arsj ■ (1) 

r,s,s' 

Here we use gauge in which tunneling amplitudes Qi and 
flrss' are real. As we noted above the key point of our 
work is that we consider the tunneling transitions accom- 
panied by spin flips . These are generated by the second 
term in Eq. (|l|) due to g-factor difference between the 
dot and the right contact. The mechanism generating 
such transitions is similar to that of spin scattering by 
nonmagnetic impurities in semiconductors Due to 
spin-orbit interaction, relatively strong in the right con- 
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tact in our case, the orbital and spin states of the electron 
in the right reservoir are mixed, resulting in effective g- 
factors for the electrons there to be different from 2. The 
eigenstates of Hn are represented by Kramers doublet, 
IV'r,s=i/2) '^r\ T)+Wr| i) and a Kramers conjugate state 
|V'r,s=-i/2}: where Ur and are functions of spatial coor- 
dinates only. We have assumed the spin orbit coupling in 
the left reservoir and the dot is much weaker (17 w 2), so 
that we can neglect by the spin-orbit mixing effect there. 
In order to evaluate the tunneling matrix elements for 
transitions from the dot to the right reservoir, given by 
the second term in Eq. ([T]), one can utilize Bardeen's 

formula §: Vlrss- = l/2m / d5 • {<i)lVipr,s' -A,s-^<P*s). 
where the integral is over any surface lying entirely within 
the tunneling barrier, separating the dot and the right 
reservoir, and the wave functions (ps (state with spin s in 
the dot, \4)s) = \4>)\s)) and ^/'r.s' are smoothly continued 
under the barrier; m is electron's mass and h = 1. It is 
obvious that the states tps under the barrier are still spin- 
orbit mixed due to the continuity condition. Therefore 
the tunneling matrix elements, corresponding to the tran- 
sitions from the resonant level to the right reservoir with- 
out spin flip, are Q.rss = l/2m J dS ■ {<p*Vur — UrVcj)*), 
and the matrix elements of transitions accompanied by 
spin flips are 0.rss — l/2m ^ dS ■ {4'*Vvr — WrV0*); 
s = —s. For relatively small deviations of g factor in 
the right reservoir from 2, \v\ ^ 0(|AgM|), /S.g — g — 2 
and so the two transition amplitudes are related as 
l^^rssl ~ 0{\Agflrss\)- For Ag > 1, the two components 
Ur and Vr are of the same order of magnitude and so 

In this work we are interested in spectral properties 
of the tunneling current and calculate its spectral den- 
sity. Typically calculations of this sort involve evalu- 
ation of the two-particle Green's functions, which is a 
quite formidable task in non-equilibrium case, beyond 
the applicability of the linear response theory. Instead, 
we adopt an alternative approach developed in Refs. , 
that allows one to evaluate the transport rate equations 
from the microscopic Hamiltonian. In this letter we show 
that one can obtain the spectral density of fluctuations 
from these equations as well (see below) . In the following 
we outline our calculation of current spectral density and 
analyze the obtained results. 

We construct the time dependent wave function of the 
system as 



\^{t)) = {&o(i) + Y^[bis{t)alais + bis{t)al 

l,s 

-Y^[birs{t)al,ais + bir-s{t)al^ais] + ... }|0), (2) 



l.r,s 



where the "ground" state |0) corresponds to the situation 
when all states below Fermi energy in the left contact are 
filled, while all states above Fermi energy in the right con- 
tact are empty. The above wave function is a superposi- 



tion of all possible particle-hole combinations that can be 
generated by the Hamiltonian H; note that H conserves 
the total number of particles in the system. Thus the first 
term in is the amplitude of the unperturbed state, i.e., 
when no excitations in the system is present, the second 
term describes a state in which a hole is created in the 
left reservoir and a particle with the same spin occupies 
the resonant level, etc. The above wave function satisfies 
the Schrodinger equation = i?!*). 

In order to describe transport in our model we intro- 
duce probabilities for the dot to be empty or occupied, 
provided that a certain number of electrons have passed 
through the junction. The level can be either empty, 
with probability a^^^, where the subscript aa indicates 
that there is no electrons in the dot and the superscript 
n describes that n electrons have arrived in the right 
reservoir /collector, or the level can be filled with prob- 
abilities (t"j and it"^, where bb indicates that the lower 
Zeeman sublevel s = 1/2 is filled, while cc stands for the 
upper Zeeman sublevel s = —1/2 being filled. Occupa- 
tion of both Zeeman levels in the dot by two electrons is 
prohibited in our model by the infinite charging energy 
U; see Refs. for detailed discussion. We also intro- 
duce the off-diagonal elements describing coherent 
superpositions of states on the upper and lower Zeeman 

levels of the electron in the dot. ct" 's are related to the 



wave function |*} as (t°^ ^ |6o|^ cr^b = Ei.s=i/2 l^'^l^ ' 

121 



'bb 

T,Ls=-l/2 = 'El,r,s[\blrs\^ + l&i,r,5|1, CtC. 

Following steps of Refs. |^,^ one derives the rate equa- 
tions for the density matrix a from the Schrodinger equa- 
tion for the wave function j'S'). These rate equations for 
a general case are presented in fl . One finds for our case: 



(^bb = -1 R(^bb + ^ Lf^aa Y~ ^ <'^' ' 



<^bc = ^^O-bc - I R'^bc - 



AT, 



AT, 



i'^bb + ^cc) 



(3a) 
(3b) 

(3c) 

(3d) 



2Trnlj^{es)pL,Ries) and ATr 



Here Tl.r 

2TTnR{es)SD,R{es)pR{£s), where we denote D,rss = ^r, 
rirss = S^R- In derivation of Eqs. (^ we assumed that 
the coupling constants il's and the densities of states 
p's are weakly dependent on energy, and so pl,r{^s) ~ 
Pl,r{£s), ^l,r{£s) = f2L,i?(es) and thus rates r^^L for 
the electrons tunneling into and out of the dot are inde- 
pendent of energy. We also assumed that Tj^ , > AF^; 
and the bias voltage condition, V ^ Tl r, which is essen- 
tial for derivation of Eqs. (^. We point out that Eqs. (|^) 
are derived in the limit of small 6Qr, therefore terms of 
order S^Ir and higher are neglected in (||). One sees from 
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Eqs. that similar to Bloch equations the two Zeeman 
levels in the dot are coupled with each other by the ofF- 
diagonal terms ( "coherences" ) due to spin-flip transitions 
through continuum with the rate AF/j. 

By summing Eqs. (^) over the number of electrons in 
the right reservoir one obtains the "standard" Bloch-type 
equations for the reduced density matrix Oij — afj 
with i,j = a, b, c. These equations, which look essentially 
identical to Eqs. (||), describe the state of the resonant 
level independently of the states of the reservoirs. 

From Eqs. one can derive the dynamics for the ex- 
pectation value of the tunneling current in the right (left) 
reservoir, (//{.lW) = *e(*(0|[i?, ^/?,L]|*(t)), where H 
is the total Hamiltonian and N^^l are the operator of 
the electron (hole) numbers is the right (left) reservoirs. 
For instance, by using Nfj — s '^Is^rs one easily 
finds from Eq. (H) that the average current in the right 
reservoir can be written as {Iu{t)) = e{Nii{t)), where 
(Nr) = E„"«a + <\ + <c)- Using Eqs. (I for a", 
one can sum over n thus obtaining 

{I nit)} = eVRlabbit) + crcc(t)] + eATnlabcit) + actit)] ■ 

It is easy to check that the transient behavior of the aver- 
age current is an oscillatory one (due to coherence terms 
ate) with frequency equal to E and approaching sta- 
tionary value 



2eTLTR (E^ 



{2Tl + Tn) [E^ + Tl) - 2Ffl AF^ 



(4) 



Actually one measures a circuit current, given by {I{t)) — 
ct{lR{t)) + P{lL{t)) where the coefficients < a,/3 < 1 
and a + P — 1 are depending on a circuit geometry . 
Obviously, the stationary current, {I{oo)) is independent 
of a and (3, Eq. (^). Yet, the transient current, {I{t)) and 
so the current spectral density are depending on a circuit 
geometry. For simplicity we consider such a case where 
only the collector current is measured (a = = 0). 
In order to evaluate current spectral density , Si{lj) — 
dr cos{uJT){I{t)I{t + r)), from rate equations (||), 
we utilize MacDonald's formula, that relates Si to the 
dispersion of charge accumulated on the collector (right 
reservoir) lioll: 



Si{lu) = e^w/TT / dtsm{ujt){N]^{t)) 



(5) 



The dispersion for the number of electrons in the 
right reservoir can be found from Eqs. (^) as {N'^) = 
J2n^'^i^2a + '^bb + ^cc)- The later must be averaged 
with respect to the stationary state distribution in the 
dot. Evaluation of the above sum is tedious but can be 
performed exactly from Eqs. (||). Here we quote the fi- 
nal expression for spectral density Si {to). The general 
result is rather cumbersome. In the region of interest, 
E > Fi,Fij > AFj^, expanding Si in powers of AFjj up 
to 0(AF|,), we obtain: 



Siii^) 



4Fi 



TT 2Fl + Tr (2Fi + TrY + ^2 



tE'^ {2Ti 



TrY Tl + {u;- EY 



(6) 



The spectrum ^ is shown in Fig. 2. The first term 
in (H) is the shot noise approaching the "Schottky" limit 
Si — e{I)/n for uj :s> F^, F^. For frequencies uj <Tl,Tr 
there is a dip in the spectrum - the result merely consis- 
tent with Refs. The second term, representing a 
distinct peak, arises due to spin-flip transitions between 
the Zeeman-split sublevels in the dot. It is roughly of 
Lorentzian shape centered approximately a,t lu = E and 
having width Tr. Importantly, the width of the peak is 
governed by F^, while the width of the dip is controlled 
by Fi for Fi ^ F^. Such condition guarantees that the 
peak is sufficiently distinct and thus can be resolved. A 
similar situation takes place in case of a current tunnel- 
ing through a double well structure where a peak 
in the fluctuation spectrum appears to be located at the 
tunneling frequency for the double well structure. 

The ratio of the peak's height to the noise pedestal 
(the signal to noise ratio) given by Eq. ^ is S/N — 
ATlATyE^TR{2TL + Tr). The 5 to ratio can be 
significantly increased in heterostructures with greater 
.91 ~ .92 difference, and thus greater spin transition rate 
ATr, or in asymmetric SETs with F^ ^ Tr. 
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Fig. 2: Power spectrum of tunneling current fluc- 
tuations; Eq. (§. Here E = 50Tr, ATr = OATr, 
Tl = 20Tr. 

The above described spin-coherence mechanism can be 
used for single spin detection. Suppose that nuclear spins 
in the dot are polarized. From Eqs. (^), (^ one can eval- 
uate the orders of magnitude for parameters needed for 
observation of the distinct peak in the fluctuations spec- 
trum in magnetic field generated by a few nuclear spins. 
The width of the peak in (^) is defined by the value of 
current through the structure, F^ ~ {I)/e for Tl > Tr. 
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Therefore in order to resolve a peak due to spin flip tran- 
sitions one needs to satisfy condition E > F^j, though E 
should not be too great as signal to noise ratio decreases 
with growth of E. Assuming that the Zeeman splitting 
E is solely due to hyperfine coupling, which is typically 
of order 10^ MHz per single nuclear spin, the measurable 
tunneling current through the structure would be of or- 
der 100 pA. This number is well within the capabilities 
of today's single-electronics. We thus conclude that the 
interference effect in resonant magneto-transport can be 
used for detection of polarization produced by < 10 nu- 
clear spins. The sensitivity of such measurements might 
be higher than that of any existing experimental setup 



We argue that the effect considered in the present work 
can explain coherent oscillations with Zeeman frequency 
in the tunneling current, which have been observed in a 
set of scanning tunneling microscope (STM) experiments 
[^5|. Impurities at semiconductor surfaces are known to 
form resonant levels that can influence the STM tunnel 
current |l^ . In contrast to previous attempts to explain 
the experiments 15 1, we suggest that it is not the impu- 
rity spin but the current itself generates coherent oscilla- 
tions due to the tunneling of electrons with spin flip via 
the resonant level formed by the impurity. 

We must note, however, that the above condition 
E > Til for the observation of a narrow peak in the spec- 
tral density requires the resonant current {Ir/s) to be 
less than the oscillation frequency E. On the first sight 
this condition seem to contradict to our statement that 
the above discussed interference mechanism can naturally 
explain experiments p5| . Indeed, this peak has been de- 
tected in the opposite regime, namely I /e > E. Yet, in 
the experiments |l5| the tunnel current flows through a 
cluster of impurities having a number of resonant levels. 
As a result the current through a single resonant level 
{Ir) represents only a fraction of the total STM current 
and therefore is considerably smaller than the total STM 
current / [0. Thus the width of the peak in the noise 
spectrum in experiments |l5| must be defined by a cur- 
rent through a single resonant level rather than by the 
total STM current. Therefore we argue that the condi- 
tion I/e > E can be met with no contradiction to the 
principal conclusions of this work. 

Finally we emphasize that the noise spectrum for the 
noise in the circuit current (/c = alji + PIl) can also be 
calculated using approach developed in this work. Using 
charge conservation, II — Ir. = Q, where Q is charge in 
the dot, one obtains a simple relation for the noise spec- 
tra of the tunnel currents through right and left junc- 
tions: Si^iuj) = aSi^iuj) + (3Si^{uj) - al3uj^SQ{Lu). The 



spectrum of charge fluctuations in the dot, Sq, can be de- 
rived from the rate equations (|^) by calculating station- 
ary charge auto-correlation function, say {(TaaiO)o'aa(t)). 
This calculation, however, will alter the results only 
quantitatively and thus will be rendered to future work 
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